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Abstract Transverse electron focusing in two-dimensional electron gases (2DEGs) 
with strong spin-orbit coupling is revisited. The transverse focusing is related to the 
transmission between two contacts at the edge of a 2DEG when a perpendicular 
magnetic field is applied. Scanning probe microscopy imaging techniques can be 
used to study the electron flow in these systems. Using numerical techniques we 
simulate the images that could be obtained in such experiments. We show that hy- 
brid edge states can be imaged and that the outgoing flux can be polarized if the 
microscope tip probe is placed in specific positions. 



1 Introduction 

During the last decade, a tremendous amount of work has been devoted to manipu- 
late and control the spin degree of freedom of the charge carriers U . It was quickly 
recognized that the spin-orbit (SO) interaction may be a useful tool to achieve this 
goal. This is due to the fact that the SO coupling links currents, spins and external 
fields. Using intrinsic material properties to control the carrier's spin would allow 
one to build spintronic devices without the complication of integrating different ma- 
terials in the same circuit JT). The challenging task of building spin devices based 
purely on semiconducting technology, requires one to inject, control and detect spin 
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polarized currents. During the last years a number of theoretical and experimental 
papers were devoted to the study of the effect of SO coupling on the electronic, mag- 
netic and magnetotransport properties of 2DEGs (see J2) and references therein). 
The nature of the SO coupling in these systems is due to the Dresselhauss and the 
Rashba mechanisms, the latter being the dominant effect in several cases Q. In ad- 
dition, the Rashba coupling has the advantage that its strength can be changed when 
a gate voltage is applied to the heterostructure, opening new alternatives for device 
design 0. 

In many transport experiments in 2DEG with a transverse magnetic field, includ- 
ing quantum Hall effect and transverse magnetic focusing, the SO coupling plays 
a central role. The transverse focusing consists basically in injecting carriers at the 
edge of a 2DEG and collecting them at a distance L from the injection point. The 
propagation from the injector / to the detector D is ballistic and the carriers can be 
focalized onto the detector by means of an external magnetic field perpendicular to 
the 2DEG. The field dependence of the focusing signal is essentially given by the 
transmission from / to D. In a semiclassical picture, the trajectories that dominate 
the focusing signal are semicircles whose radius can be tuned with the external field. 
The new scanning technologies developed in |5]|6) can be used to map these trajec- 
tories. The scanning probe imaging techniques consist in perturbing the system with 
the tip of a scanning microscope and plotting the transmission as a function of the 
tip position. The transmission change is a map of the electron flow. In this paper 
we first revisit the theory of transverse electron focusing in systems with strong 
SO coupling and interpret the results in terms of a simple semiclassical picture Q. 
Then, we use numerical techniques to simulate the images that could be obtained 
with scanning probe microscopy experiments. We show that hybrid edge states can 
be visualized and that the outgoing flux can be polarized if the microscope tip probe 
is placed in specific positions 



2 Transverse electron focusing in presence of strong spin-orbit 
coupling 

The Hamiltonian of a 2DEG with Rashba spin-orbit coupling is given by 

H = i ( P ? +P ^ + f Vy°x- p * a y)- \g^BO z B z +V{x) (1) 

here P ri =p-q — (e/c)A^ with and being the ^-component of the momentum 
and vector potential respectively, a is the Rashba coupling parameter, g is the ef- 
fective g-factor, {cijj} are the Pauli matrices and V(x) describes the potential at the 
edge of the sample. In what follows, we use a hard wall potential: V(x) =0 forx > 
and infinite otherwise. For convenience we choose the vector potential in the Landau 
gauge A=(0,*B z ,0). 
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Far from the sample edge (x 0) the eigenvalues and eigenfunctions of Hamil- 
tonian (Q]i are well known (SJ. The SO coupling breaks the spin degeneracy of the 
Landau levels. The spectrum is given by 
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E±=h(D c nT\jEl+^jj 2n, (2) 

where n > 1, co c = e\B \ /m*c is the cyclotron frequency, l c = (Ti/ma) c )? is the mag- 
netic length, and Eo = ha> c /2 — g^igB-Jl is the energy of the ground multiplet corre- 
sponding to n=0. The eigenfunctions for n > 1, written as spinors in the z-direction, 
are J9) 

and 

where L, is the length of the sample in the y-direction, (j) n (x — xq) is the harmonic 
oscillator wavefunction centered at the coordinate xq = l^k , A n = 1 + D\ and 



D n = W '^ 2n (5) 



E + V £ o + («/^) 2n 

The wave functions of the first Landau level are given by (x,y) with « = 0. 
These eigenstates have a cyclotron radius given by 

r c 2 =2(^ ± |(x-x ) 2 |^ ± ) ! (6) 

that for large n gives r 2 . ~ 2n(h/m* (O c ). We see from Eq.© that states with differ- 
ent «, and consequently different cyclotron radius, coexist within the same energy 
window. Additionally, in the limit of strong Rashba coupling or large n, D„ ~ 1 and 
the spin lies in the plane of the 2DEG. 

Equivalent results are found in a semiclassical treatment of the problem lfTTl[Tol . 
In this approach, the spin is described by a vector ifTTl S = h/2(n\{i) ,ri2(t) ,n^(t)) 
and the classical orbits are given by 

q = r±(coso±r,sina)±f) 
h 

S = sign(B,)-(=Fcoso)±f,=Fsina)±f,0), (7) 
here q is the coordinate measured from the centre of the circular orbit of radius 

If a \ 2 2E~ a 
r± = \ (i—) +-— r±r-, (8) 
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and the corresponding cyclotron frequencies are 

(0± = sign(fl z )(cj c T a/Kr±) ■ (9) 

In agreement with the quantum results obtained for large n, the spin is found to 
be in-plane pointing outwards for the smaller orbit and inwards for the bigger one 
when a positive perpendicular magnetic field B z is applied. Moreover, the Born- 
Sommerfeld quantization 11121 of these periodic orbits reproduces the exact quantum 
results of Eq.© for large «. 

The calculation of the exact edge states with the hard wall potential requires a 
numerical approach. We have shown that the semiclassical approximation can be 
extended to describe edge states in which electrons bounce at the sample edge iflOl . 
Due to the continuity of the wave function and the spin conservation at the edge, 
the two orbits with radii r + and r_ are mixed as schematically shown in Fig.l. The 
agreement between the Born-Sommerfeld quantization of the semiclassical edge 
states and the quantum results is excellent for states composed of semicircles cen- 
tered in the edge (normal incidence). In what follows, we use these semiclassical 
orbits to interpret the numerical results for transverse focusing experiments. 

The transverse focusing experiments collect electrons or holes coming from a 
point contact |[T3l[T4l into another point contact acting as a voltage probe. The car- 
riers are focused onto the collector by the action of an external magnetic field as 
schematically shown in Fig. 1. The signal measured in transverse focusing experi- 
ments is related to the transmission T between the two point contacts located at a 
distance L from each other (see Fig. 1). Typical experimental setups also include 
two ohmic contacts at the bulk of the 2DEG which are used to inject currents and 
measure voltages. The details of different configurations with four contacts have 
been analyzed in fl5l . The main features of the magnetic field dependence of the 
focusing peaks are contained in T lfl6l . Consequently, from hereon we will refer to 
the focusing signal or to T indistinctly. In the zero temperature limit we only need 
to evaluate T at the Fermi energy Ep. For the numerical calculation of T the system 
was discretized using a tight-binding model in which the leads or contacts are easily 
attached. In this approach the Hamiltonian is given by H = Hq + H$o with 

ffo = ^ £ ac\a c na ~ H ^nmc\ a c ma +h.c. . (10) 

«,(7 <«,m>,c 

Here cj, a creates an electron at site n with spin a (f or j in the z direction) and 
energy £,j =4?— (7g/iB,B z /2, t = — Tr /2ot*Oq and qq is the lattice parameter which is 
always chosen small compared to the Fermi wavelength. The summation is made on 
a square lattice, the coordinate of site n is n x x + n y y where x and y are unit vectors 
in the x and y directions, respectively. The hard-wall potential V(x) is introduced 
by taking n x > 0. The hopping matrix element t nm connects nearest neighbors only 
and includes the effect of the diamagnetic coupling through the Peierls substitution 
ifTTl . For the choice of the Landau gauge f„(„+y) = t exp (m x 2nij) /0q) and t„/ n+ -^ =t, 
=OqB z is the magnetic flux per plaquete and <j>o=hc/e is the flux quantum. The 
second term of the Hamiltonian describes the spin-orbit coupling, 
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Fig. 1 Panel(a) Transmission coefficient between the contacts / and D as a function of the applied 
perpendicular magnetic field in the presence of strong Rashba spin-orbit coupling (qualitative). 
Relevant semiclassical orbits for three different focusing conditions are shown in panels (b),(c) 
and (d). 



H SO - {^ c lf{n+x)l ~ K c l[ c (n+x^ + (11) 



e in x 2K<j>/<l) 



h c l] c (n+y)[-K c l[ c (n+y)] } + h 



where X x ~ a/2ao and Xy = — ia/2«o- In what follows we use the following val- 
ues for the microscopic parameters: ciq — 5nm, m* = 0.055rao — here mo is the free 
electron mass — and Ep = 23meV. These parameters correspond to InAs based het- 
erostructures with a moderate doping. We use different values of the SO coupling 
parameter a as indicated in each case. 

The two lateral contacts, / (injector) and D (detector) are attached to the semi- 
infinite 2DEG described by Hamiltonian ( fTTT i. Each contact is an ideal (with a = 0) 
narrow stripe of width Nquq. They represent point contacts gated to have a single 
active channel with a conductance 2e 2 /h, for details see Q. 

To obtain the transmission between the two contacts we calculate the Green func- 
tions between the sites of the injector and the sites of the detector. As the spin 
is not conserved, the Green function between two sites i and j has four compo- 
nents tyo.jo'- Fi rst tne propagators of the system without the contacts are obtained 
by Fourier transforming in the y-direction and generating a continuous fraction for 
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each k. Having these propagators, the self energies due to the contacts can be easily 
included using the Dyson equation [17]. The transmission is then obtained as 



~h 



r (2)yR r (l) 



(12) 

=E F 



here Sf R and Sf A are the retarded and advanced Green function matrices with ele- 
ments - a i and 'S^ - c , . The matrices are given by the self-energy due to 

contact /, — ]\Lf— £f] where Ef wdZf are the self-energies matrices of the re- 
tarded and advanced propagators respectively. Note that the definition of T includes 
the spin index. 

A typical T vs. B z signal for strong spin-orbit coupling is shown in Fig[TJa). 
A splitting of the first focusing peak is clearly observed \J\. Notably, there is no 
splitting of the second peak. These results can be easily interpreted in terms of the 
semiclassical picture given above. From all the semiclassical orbits that connect the / 
and D contacts, the ones that give the largest contribution to T are those with 2r± = L 
17] H31 . When the applied magnetic field B z is increased the cyclotron radii are 
reduced as Br 1 and the first maximum in the transmission is found when r_ (B z ) = 
L/2 as schematically shown in Fig[T](b). There Oi is the electron path between 
/ and D, this path is a semicircle of radius r_. For this field, indicated as B- = 
B\ i, the electrons that flow out of the injector in the O2 orbit do not arrive to the 
detector since r + (B\ 1) > L/2. Furthermore, the two orbits 0\ and O2 correspond 
to electrons injected with spin down or up in the y-direction, respectively. Note 
that due to the SO coupling, the spin rotates along the orbit. It is convenient to 
split the total transmission in the four contributions T a p corresponding to electrons 
injected with spin a and collected with spin j3. The total transmission can be put 
as T = T m + T lu ( + T<j u + T^d and for B- = B 1 , 1 the total transmittance is dominated 
by the contribution 7j„. When B z is increased over Si i, r_(5i,i) < L/2 and T 
decreases. The next maximum is reached for B Z =B\2 when r + {B\o)=L/2 and the 
relevant orbit is O2 as shown in FigQ](c). For this focusing field the transmission is 
dominated by T m i- 

The next maximum in T is found when 5- = 52 an d corresponds to the situation 
shown in FigQ](d). This focusing condition is due to the semiclassical trajectories 
with one intermediate bounce at the edge of the sample. In this case the two possible 
paths 0\ and O2 contribute to the transmission. Electrons leaving the injector with a 
given spin arrive at the detector with the same spin projection. Accordingly, the total 
coefficient T is dominated by T uu ^-T c [d- Clearly, B2 is the magnetic field for which 
2(r_+r + ) =L holds. In agreement with the exact numerical result, by extrapolating 
the semiclassical picture shown in FigJT] one finds that the peaks that are split due 
to Rashba interaction are those in which the number of bounces is even (or zero). 
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Fig. 2 (a) Total focusing transmission coefficient T versus applied perpendicular magnetic field 
B z . (b) Spin resolved transmission coefficients versus B : . We used Ef = 23meV, m* =0.055mo, 
a =7meVnm, j8 =0, L= l.Sflm and the width of the contacts is IQnm. (c) Schematic of an SPM 
imaging procedure. 



3 Imaging Techniques in Transverse Focusing with spin-orbit 
coupling 

Scanning probe microscopy (SPM) techniques have been recently used for imaging 
the electron flow in a variety of 2DEG ballistic systems J5] |6). With this technique, 
the negatively charged tip of a scanning microscope is positioned above the 2DEG 
as schematically shown in Fig|2jc). The tip position can be changed to sweep a 
given area of the explored 2D device. The electrons under the tip are repelled and 
consequently a zone of lower electron density (or divot) is formed under the tip. 
In the simplest case the transmission (and then the conductance) between two con- 
tacts of the device is measured as the tip position changes. If the tip is located in a 
region that affects the electron path between the contacts the conductance changes 
providing a map of the electron flow in the device. The resolution of these images is 
smaller than the divot size, ||5][6) making this technique a powerful tool for studying 
nano-scale ballistic systems. 

Here, we propose the use of this technique to explore the transverse focusing in 
the presence of spin-orbit interaction lfl"8l . We simulate the effect of the tip potential 
by perturbing (increasing) the site energies £,,<j in an area of the order of 100«m 2 
centered at the tip position. The Dyson equation is used to introduce the perturbation 
and the exact propagators between the contacts / and D are calculated for each 
position of the tip. 

Figure [3] (a) shows T vs the tip's position when the perpendicular magnetic field 
is fixed to obtain the first maximum {B Z =B\ \) for a SO coupling a = 7meVnm. 
The semicircular electron path is clearly observed. In this case the T map is dom- 
inated by a drop in T along the 0\ path. A similar pattern is found for the second 
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Fig. 3 Total transmission coefficient T between the contacts 1 and 2 as a function of the probe 
position for (a) B- = B lA , (b) B- = B L2 and (c) B z = B Lr = (B lA + B L2 )/2. We used a = 7meVnm 
and the parameters given in Fig[2] 



transmission maximum {B Z =B\2) as shown in Fig(3j(b). In this case, the drop in T 
is due to the scattering induced by the tip of electrons travelling along the O2 path. 
A slightly different situation is found when B z is fixed in between B\ \ and B\ 2 as 
shown in Fig[3](c); although the variation is also dominated by a drop (dark area), 
T increases at the two sides of the minimum. The observation of these two lobes 
shows that the tip, when placed at those positions, modifies the electron flow mak- 
ing a non-focalized electron path — 0\ or O2 in Fig[T](b)-(c) — to contribute to the 
transmission. 

More interesting are the imaging results obtained when the field is fixed at the 
second focusing condition: B Z =B2- As mentioned above, in this case T is domi- 
nated by the electron's orbits with one bounce at the sample's edge. For this field 
the largest contributions to the transmission coefficient are T uu and T c / C / , and the cor- 
responding focusing peak is unsplit. In Fig|4]and Fig[5]the results for this case are 
shown for a = 7meVnm and 15meVnm, respectively. Panel (a) shows T m as a func- 
tion of the position of the microscope probe. The change in the transmission in this 
case clearly shows that the electrons injected with spin up (in the y-direction) leave 
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Fig. 4 We plot (a) T uu , (b) T c u, (c) T and (d) P as a function of the probe position for B z =Bi. We 
used a = 7meVnm and the parameters given in Fig|2] 



the injector in the bigger orbit, rebound and then arrive to the detector in the smaller 
orbit with spin up — see O2 in Fig[TJd). In panel (b) the transmittance T c id is shown 
— see 0\ in Fig[TJd) - and in panel (c) the total transmission coefficient is presented. 
As T uc i and T c / U are very small, the total transmission is essentially given by the sum 
of the contributions shown in panels (a) and (b). Experimentally these two contri- 
butions could be distinguished by selecting the spin of the injected carriers. In fact, 
a combination of an external in-plane magnetic field in the y-direction and an ap- 
propriate gate voltage in the point contacts can be used to filter spins in the injector 
or detector lfl3ll . Selecting the spin of the injected electrons would make it possible 
to separate the two trajectories — (a) and (b) in Figs. [4] and [5] — and obtain a direct 
visualization of the two orbits split by the spin-orbit coupling. Conversely, selecting 
the spin in the detector D, the transmissions T + = T uu + T c [ u and T = T m / + of 
carriers arriving at D with spin up and down, respectively, could be measured. In 
terms of these quantities, we define the polarization P of the transmitted particles as 
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Panel (c) and (d) of Figs. |4] and [5] show the total transmission coefficient T and 
the polarization P as a function of the tip position. The two semicircular electron 
paths including the rebound at the edge are visualized in the T map. In our sim- 
ulations the smaller and the bigger electron paths are not easily resolved in total 
transmission coefficient map except for the largest SO coupling case and for the tip 
close to the bounce position — see Fig|3c). There, an appreciable fall (about 50%) 
of the transmission in the two rebound positions indicates that, when the probe is 
positioned there, the contribution to T of one of the two possible electron paths {0\ 
or O2) is being suppressed. If Oi is being suppressed, the electrons arriving to the 
detector will have spin up. On the other hand, if O2 is suppressed only spin down 
electrons will arrive to the detector. This means that one can select the spin polar- 
ization of the outgoing carrier flux by changing the tip position a few nanometers 



Imaging transverse electron focusing in 



11 



as shown in FigHfd) and Fig(3d). Notably, the effect is also clearly observed in the 
case of the smaller SO coupling despite of the fact that the total transmittance T 
does not resolve the two orbits. 



4 Summary and Conclusions 

We have discussed a microscopy imaging technique for the case transverse electron 
focusing in 2DEGs with strong Rashba coupling. The main results can be summa- 
rized as follows: 

;.- The existence of two different cyclotron radii splits the first focusing peak 
onto two sub-peaks, each one corresponds to electrons arriving to the detector with 
different spin polarization along the direction parallel to the sample's edge. 

ii- The images of the electron flow for focusing fields corresponding to the first 
two sub-peaks are very similar and consequently, for this case, the technique can 
not clearly distinguish the two type of orbits. 

Hi.- When the external field is fixed between the focusing fields of the two sub- 
peaks, B z = (B\ i +B\ t z)/2, the transmission map shows a structure that indicates 
the presence of the two orbits. 

iv. - For the second focusing condition, and for the case of strong Rashba cou- 
pling, the technique can resolve the two orbits when the microscope tip is placed 
close to the rebound position. 

v. - For the case described in the previous point, the microscope tip can be used to 
polarize the electron flux arriving at the detector. The direction of the polarization 
can be reversed by changing the tip position a few nanometers. 

Finally, we would like to emphasize a few points: a.- Interference fringes, charac- 
teristic of the quantum ballistic transport regime, are observed in all the T maps; b.- 
For the properties studied here, replacing the hard wall potential V(x) by a more re- 
alistic parabolic potential does not change the main properties of the system 111 911201 . 
Therefore, our results should correctly describe the images that could be obtained 
in heterostructures defined by gates, c- The competition between the Rashba and 
the Dresselhauss couplings leads to interesting features in the focusing signal and 
needs to be considered for interpreting imaging results in systems where these two 
SO interactions are present. These results will be presented elsewhere. 
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